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ABSTRACT 

We present a study of the X-ray properties of a sample of six nearby late-type spiral galax- 
ies based on XMM-Newton observations. Since our primary focus is on the linkage between 
X-ray emission and star formation in extended, extranuclear galactic disks, we have selected 
galaxies with near face-on aspect and sufficient angular extent so as to be readily amenable to 
investigation with the moderate spatial resolution afforded by XMM-Newton. After excluding 
regions in each galaxy dominated by bright point sources, we study both the morphology and 
spectral properties of the residual X-ray emission, comprised of both diffuse emission and 
the integrated signal of the fainter discrete source populations. The soft X-ray morphology 
generally traces the inner spiral arms and shows a strong correlation with the distribution of 
UV light, indicative of a close connection between the X-ray emission and recent star for- 
mation. The soft (0.3-2 keV) X-ray luminosity to star formation rate (SFR) ratio varies from 
1 — 5 x 10 39 erg s _1 ( M Q yr -1 ) -1 , with an indication that the lower range of this ratio relates 
to regions of lower SFR density. The X-ray spectra are well matched by a two-temperature 
thermal model with derived temperatures of typically ~ 0.2 keV and ~ 0.65 keV, in line with 
published results for other normal and star-forming galaxies. The hot component contributes a 
higher fraction of the soft luminosity in the galaxies with highest X-ray/SFR ratio, suggesting 
a link between plasma temperature and X-ray production efficiency. The physical properties 
of the gas present in the galactic disks are consistent with a clumpy thin-disk distribution, 
presumably composed of diffuse structures such as superbubbles together with the integrated 
emission of unresolved discrete sources including young supernova remnants. 
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1 INTRODUCTION 

The first detailed studies of the X-ray emission from nearby late- 
type galaxies were carried out with the Einstein observatory and 
revealed both individual luminous X-ray sourc es, as well as a n 
extended, often complex, underlying structure jFabbiand[l98^) . 
In a number of nea r edge-on star-forming systems, such as M82 
jWatsonetai]|l984h . evidence was also found for hot outflowing 
winds. The superior spatial resolution and enhanced soft X-ray re- 
sponse of ROSAT greatly extended this work, leading to the de- 
tection of substantial numbers of discrete sources as well as es- 
tablishing a much clearer p icture of the diffuse emission compo- 
nents (e.g. iRead et al]| 19971) . More recently the advent of XMM- 
Newton and Chandra has extended our knowledge of the discrete 
X-ray source populations and the ultrahot ISM in nearby galax- 
ies to u nprecedented luminosity and surface brigh tness thresh- 
olds (e.g. IStrickland et alJl2004atlColbert et al.ll2004 iPietsch et al.l 
2004lFabbianoll2006l : lYang et afcollBogdan & Gilfanovll2008T 
Bauer et al.ll2008h . 



The diffuse X-ray emission associated with the disks of spi- 



ral galaxies is believed to arise in gas heated by shocks generated 
as a result of supernova explosions and stellar w i nd interactions 
dChevalier & Clegd 1 19851 : IStrickland etafl l2000l : ISuchkov et ail 
ll994lGrimeset alj|2005r) . This is supported by analysis of the X- 
ray spectra of individual shell-type supernova remnants (SNRs) of 
reasonable age, which are often characterised by plasma tempera- 
tures in the range 0.2-0.8 keV, similar to the value s obtained when 
studyin g larger-scale structures in galactic disks. Strickland et al. 
( I2004ah showed that the luminosity of the diffuse X-ray emis- 
sion observed in star-forming galaxies is proportional to the rate 
of mechanical energy injection into the I SM from young stars. 
iMas-Hesse, Oti-Floranes & Cervine] |2008) further found that the 
ratio of the soft X-ray luminosity to far infrared luminosity, taking 
the latter as a proxy for the star-formation rate (SFR), is strongly 
dependent on the efficiency with which mechanical energy is con- 
verted into soft X-rays and also the evolutionary status of the star 
formation episode. This all argues for a close connection between 
soft X-ray emission and star formation in late-type galaxies. 

A recent study of a sample of 12 nearby spiral galax- 
ies using Chandra dTvler et all 120041) has shown strong cor- 



2 R.A.Owen & R.S. Warwick 

Table 1. Properties of the sample of galaxies. 



Galaxy 


Hubble Type a 


Distance b 
(Mpc) 


D 25 a 
(') 


Inclination b 

(°) 


Foreground Nh c 
(10 20 cm- 2 ) 


Nucleus 


NGC300 


SA(s)d 


2.0 


21.8 


30 


3.6 




M74 


SA(s)c 


11.0 


10.5 


7 


4.8 




NGC3184 


SAB(rs)cd 


11.6 


7.2 


< 24 


1.0 




M51 


SA(s)bc 


8.4 


7.3 


20 


1.8 


Sy 2 


M83 


SAB(s)c 


4.5 


12.9 


24 


3.9 


Starburst 


M101 


SAB(rs)cd 


7.2 


26.9 


18 


1.1 





a - Hubble type and the major-axis diameter (D25) from the RC3 catalogue Ide Vaucouleurs et al J 199 lh . 

b - See j|4.2| for references. 

c - Njj values based on iDickev & Lockmanl fl99(il) . 



Table 2. Details of the XMM-Newton observations of the five new galaxies in our sample. 



Galaxy 


Observation ID 


Start Date 


Filter a 


Target co-ordinates b 




Useful exposure (ks) 






(yyyy-mm-dd) 


pn/MOSl/MOS2 


RA (J2000) 


Dec (J2000) 


pn 


MOS 1+2 


NGC300 


0112800201 


2000-12-26 


M/M/M 


00 ,l 54 m 55.0 s 


-37° 41' 


00" 


29.4 


65.1 




0112800101 


2001-01-01 


M/M/M 








40.4 


88.4 




0305860401 


2005-05-22 


M/M/M 








27.2 


64.9 




0305860301 


2005-11-25 


M/M/M 








30.4 


68.2 


M74 


0154350101 


2002-02-02 


T/T/T 


0l' l 36 m 23.9 s 


+15° 45' 


13" 


23.3 


49.9 




0154350201 


2003-01-07 


T/T/T 








23.5 


50.2 


NGC3184 


0028740301 


2001-11-02 


T/T/T 


10 h lg m 40 qs 


+41° 25' 


11" 


19.9 


49.9 


M51 


0112840201 


2003-01-15 


T/T/T 


13 ft 29 m 51.9" 


+47° 10' 


32" 


18.8 


40.6 




0212480801 


2005-07-01 


M/M/M 








20.4 


56.7 




0303420101 


2006-05-20 


T/T/T 








32.3 


77.3 




0303420201 


2006-05-24 


T/T/T 








14.1 


43.0 


M83 


0110910201 


2003-02-18 


T/M/M 


13' l 37 m 00.4 s 


-29° 52' 


04" 


18.3 


38.2 



a - T = thin filter, M = medium filter. 

b - Assumed position of the galactic nucleus. 



relation between soft diffuse X-ray emission and sites of re- 
cent star formation in spiral arms traced by mid-infrared and 
Ha emission. A similar investigation has not been completed us- 
ing XMM-Newton, but a nu mber of studies of i ndividual spiral 
galax i es have been reported { Pietsch et alj|200ll; iTakahashi etall 
l2004lTrudolvubov et al1l2005l : IWarwick et al.ll2007h . which again 
demonstrate the strong linkage between tracers of star formation 
and diffuse X-ray emission. Furthermore, a dual-component ther- 
mal model, with characteristic temperatures of ~ 0.2 keV and 0.6- 
0.7 keV, very often provides a good description of the spectrum 
of the diffuse X-ray emission. Studies of the point source X-ray 
luminosity function (XLF) in spiral g alaxies dTennant et al.ll200ll ; 
ISoria & Wull2003l : ICoibert et alfeOoA have shown that the bright- 
est point sources contain the bulk of the integrated point source 
luminosity, suggesting that, provided a sufficient number of the 
brightest point sources are excluded, it should in principle be possi- 
ble to probe the underlying structure of the galaxy to relatively low 
levels of surface brightness. 

In this paper, we focus on the spectral and morphological 
properties of the diffuse X-ray emission emanating from the ex- 
tended disks of late-type normal galaxies, as deduced from XMM- 
Newton observations. In Sj2] we identify a sample of nearby galax- 
ies with close to face-on disks of sufficient angular extent so as 
to be readily amenable for study with XMM-Newton. We then out- 
line the data reduction methods employed (0. Using image anal- 
ysis, we sub-divide the total X-ray luminosity of each galaxy into 
the contribution of spatially resolved bright sources and a residual 
component comprising the integrated emission of unresolved point 



sources and the diffuse galactic emission (Sj4](. On the basis of pub- 
lished source luminosity functions, we estimate the likely contri- 
bution of truly diffuse emission to this residual signal. We go on 
to investigate spectral properties of both the bright source and the 
unresolved (residual) components ( Ej5j. We also consider the mor- 
phology of the residual emission in the context of the star formation 
evident in the galactic disks (ij6}. Finally we discuss the implication 
of our results (SjTJ and briefly summarise our conclusions (Sj8}. 



2 GALAXY SAMPLE AND OBSERVATIONS 

In this paper we investigate the large-scale spatial and spectral 
properties of the X-ray emission emanating from the extended 
disks of six nearby late-type spiral galaxies. The galaxies in 
question are NGC300, M74, NGC3814, M51, M83 and M101, 
brief details of whic h are given in Table Q] In a previous paper 
dWarwicket"aiT2 007; hereafter W07), we reported the results of an 
XMM-Newton study of one of these galaxies, namely M101, and 
here we apply a similar methodology to extend the sample, thereby 
allowing the intercomparison of the galaxy X-ray characteristics 
based on XMM-Newton observations. 

When dealing with X-ray spectral imaging data of only mod- 
erate spatial resolution (as afforded by XMM-Newton), the pres- 
ence of a number of very bright point sources has the potential 
to obscure the properties of the underlying galaxy. Since our goal 
is to study the latter rather than the former, we have limited our 
galaxy sample to nearby face-on systems of Hubble type Sbc-Sd 
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with major-axis D25 extent greater than 7' , thereby ameliorating 
the source confusion problem. (We have excluded M33 from our 
study since this galaxy has a D25 diameter of 70' and thus ex- 
tends well beyond the EPIC field of view; M33 is the subject of 
a specific XMM-Newton programme utili sing multiple pointings 
dPietsch et al.ll2004lMisanovic et alj|2006h . 

Details of the XMM-Newton EPIC observations of the five new 
additions to our galaxy sample are summarised in Table [2] 



3 DATA REDUCTION 

The data reduction was based on SAS v8.0. The datasets were ini- 
tially screened for periods of high background by accumulating 
full-field 10-15 keV light curves. MOS data were excluded during 
periods when the 10-15 keV count rate averaged over 100 s bins 
exceeded 0.2 ct s -1 , and the pn data were similarly excised when 
the pn count rate exceeded 2 ct s -1 . As summarised in Table[2] the 
exposure times after event filtering within individual observations 
ranged from 14 to 40 ks in the pn camera (with typically slightly 
greater exposure in the MOS cameras). 

In W07, we discussed how the use of an appropriate spatial 
mask greatly reduced the "contamination" arising from the pres- 
ence of relatively small numbers of luminous point sources in the 
disk of M 101. Here we employ a refinement of this masking tech- 
nique in our analysis of the XMM-Newton observations of five fur- 
ther target galaxies. 



3.1 Image Construction 

The starting point of our analysis, for each galaxy and for each in- 
dividual observation, was the creation of pn, MOS 1 and MOS 2 
images and associated exposure maps in three energy bands: soft 
(0.3-1 keV), medium (1-2 keV) and hard (2-6 keV). In all cases, 
the pixel size was set at 4" x 4". For both the image analysis and 
the subsequent spectral analysis, we utilised single and double pixel 
events in the pn camera (pattern 0^) and single to quadruple events 
(pattern 0-12) in the MOS cameras. The resulting images were then 
flat-fielded by subtraction of a constant background particle rate 
(estimated from count rates recorded in the corners of the CCD de- 
tectors not exposed to the sky) and division by the relevant exposure 
map. At the same time, spurious data from bad pixels, hot columns 
and chip gaps were excluded. Where there was more than one ob- 
servation, the images pertaining to a given instrument and band 
were at this point combined into a single mosaiced image. For each 
band, the resultant flat-fielded MOS 1 and MOS 2 images were 
combined to give a "MOS" image and finally a "pn+MOS" image 
was also produced by summing the individual pn and mos images 
(with appropriate scaling in regions affected by chip gaps, where 
there was incomplete coverage in one of the instrument channels). 

Figures Q] & [2] show the soft-band images derived for each 
of the five new galaxies, together with the corresponding data for 
M101 from W07. In the case of M74, NGC3184, M51 and M83 our 
analysis pertains to circular regions centred on the galactic nucleus 
of diameter equal to the major-axis D25 dimension, whereas for 
NGC 300, given the relatively low X-ray surface brightness, we re- 
strict our attention to a central 10' diameter region (0.46 of the D25 
extent). Similarly for M101 the X-ray analysis reported in W07 ap- 
plies to a central 20' diameter region (0.74 of the D25 dimension). 
Hereafter we refer to these galaxy-centred circular regions as the 
"X-ray extraction regions". 



3.2 Spatial masking of bright sources 

In order to produce a spatial mask to suppress the bright source con- 
tamination it was necessary to construct a catalogue of the bright 
sources associated with each galaxy. The starting point was a pre- 
liminary source-list for each galaxy derived from the second X MM- 
Newt on Serendipitous Source Catalogue (2XMM; IWatson et al.l 
2008), encompassing those sources located within the galaxy ex- 
traction region (as defined above). Next we determined the net 
count-rate of each source in the broad-band pn+MOS imag^J 
Those sources for which the measured count rates exceeded a de- 
fined count-rate threshold (specific to each galaxy) were then in- 
cluded in the final "master" source-list for the galaxy. In the case of 
both M83 and NGC 300 two additional sources, not included in the 
2XMM catalogue, but clearly present in the images at count-rates 
above the threshold were added by hand. Similarly in the case of 
M51 one additional source was added. The number of sources in- 
cluded in each master source-list is given in Table [3] together with 
the threshold setting (quoted as the equivalent unabsorbed source 
luminosity in the 0.3-6 keV band on the basis of a source spec- 
trum consisting of power-law continuum with V = 1.7 absorbed by 
the Galactic foreground Nh). At this stage we also determined the 
count rates for each source in the soft, medium and hard bands by 
applying the above process to the appropriate pn+MOS sub-band 
images. 

The next step involved the creation of the spatial mask. To this 
end, we produced a simulated image of the point sources present 
in the broad-band pn+MOS observation. The MOS- and pn-camera 
PSFs were modelled for each energy band as a set of sub-images 
representing the PSF out to 1' from the source position. A simulated 
sub-image of each source (in three bands and two instruments) was 
then produced by scaling the appropriate PSF image so as to match 
the actual image data (in terms of the count rate contained within 
the nominal source cell). The resultant set of sub-images were then 
co-added with suitable spatial offsetting so as to simulate the sur- 
face brightness distribution in the broad-band pn+MOS image de- 
riving from bright point sources within the galaxy fielc0. A spatial 
mask was then derived from this image by applying an appropriate 
surface brightness cut. 

Figure [3] shows the simulated "bright source image" derived 
using the above procedure for M83. In this example, the bright 
source mask was obtained by applying a surface brightness cut to 
the image at a level of 0.07 pn+MOSl+MOS2 ct ks -1 pixel -1 , 
leading to a "spillover" fraction, (i.e., the fraction of the bright 
source signal not contained within the masked region) of 9%. Simi- 
lar surface brightness cuts were applied to the other galaxies result- 
ing in spill-over fractions in the range 9-12%. 

The derived spatial masks were used to separate the "bright 
source region" from the "residual emission region" in both the spa- 
tial and spectral analysis. In order to study the spatial distribution of 
the residual emission (see Sj6j, a simulated image representative of 



1 The pn+MOS 1+MOS2 count-rate above the local background was deter- 
mined within a cell of 16" radius centred on the source position. A correc- 
tion factor was then applied to allow for the signal contained in the wings of 
the point spread function (PSF) extending beyond the nominal source cell. 

2 As a final iteration, count rates were extracted from the simulated images 
within a 16" cell centred on each source position and compared to those 
obtained from the actual data, with any differences attributable to the PSF 
spreading of signal from adjacent sources. Based on this analysis, the nomi- 
nal source count rates were adjusted to compensate for this source confusion 
effect and a new version of the simulated images was produced. 
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Figure 1. Comparison of the soft X-ray and far-ultraviolet images of three galaxies. Left-hand panels: Adaptively smoothed versions of the XMM- 
Newton pn+MOS images in the soft (0.3-1 keV) band. The amplitude scaling in all cases is logarithmic with the contour levels representing factor two 
steps in the soft X-ray surface brightness. Right-hand panels: The GALEX FUV (\ e f f = 2267 A) image on the same spatial scale as the X-ray data (except 
for NGC 3184 where the UV data are from the XMM-Newton Optical Monitor in the UVW1 (ps 2680 A) band. The amplitude scaling is again logarithmic. 
The contours show the soft X-ray morphology for comparison purposes. Top panels: NGC 300. The circle has a diameter of 10'. Middle panels: M74. The 
circle has a diameter of 10.5'. Bottom panels: NGC3184. The circle has a diameter of 7.2'. 



X-ray emission from spiral galaxies 5 




Figure 2. As for Fig.[T] Top panels: M51. The circle has a diameter of 7.3'. Middle panels: M83. The circle has a diameter of 12.9'. Bottom panels: M101. 
The circle has a diameter of 20' . 
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Figure 3. Left panel: The modelled "bright source" image constructed for M83 from a set of count-rate scaled PSF sub-images. A source mask was produced 
by thresholding this image at a surface brightness level such that 91% of the source signal was contained within the resulting mask. Right panel: The soft-band 
image of the residual emission in M83 obtained from the full image by subtracting the bright source model and then applying the spatial mask. Much of the 
"contamination" due to bright sources is suppressed by this process. In both panels the circle has a diameter of 12.9'. 



Table 3. Characteristics of the bright X-ray source population in each galaxy. 



Galaxy 


X-ray Region™ 
(') 


Threshold Lx 6 
(10 37 ergs" 1 ) 


Number in 
Source List 


Number of high Lx sources 
(L x > 5 X 10 38 ergs" 1 ) 


NGC300 


10.0 


0.2 


22 





M74 


10.5 


5.0 


20 


4 


NGC3184 


7.2 


6.0 


18 


4 


M51 


7.3 


4.0 


26 


6 


M83 


12.9 


2.0 


40 


2 



a - Diameter of the "X-ray extraction region". 
b - Nominal Lx threshold applied in the 0.3-6 keV band in defining sample. 



the soft band data was subtracted from the corresponding pn+MOS 
image and the source mask imposed - see Figure [3] This process 
serves to to suppress much of the "contamination" arising from the 
brightest sources in the galaxy, including the bulk of the spillover 
into the residual galaxy region. Note that a slightly different ap- 
proach is necessary in the spectral analysis; here we take account 
of the bright source signal spilling into the residual emission region 
by including an appropriate fraction of the bright source spectral 
model in the spectral fit of the residual emission - see S]5.2I 



3.3 Spectral Extraction 

The soft-band images shown in Fig. Q]& [2] demonstrate the exis- 
tence of an extended emission component in addition to the pop- 
ulation of bright point sources. On the basis of the approach de- 
scribed earlier, we extracted the integrated spectrum of both the 
bright source region (bounded by the spatial mask) and the residual 
emission region (corresponding to the full X-ray extraction region 
less the masked area). 

The relatively large extent and complicated shapes of the 
bright source and residual emission regions makes the evaluation 
of the background subtraction somewhat more complicated than 
with most XMM-Newton data. In this context, the use of "blank- 
sky" XMM fields to produce a background spectrum would be far 



from ideal, because the sky X-ray background in the soft band 
varies from field to field, and this is a vital component of the to- 
tal background signal we wish to subtract. Instead, we extracted 
spectra from both an annulus surrounding the defined galaxy re- 
gion and also from the corner regions of each detector not exposed 
to the sky. By using an appropriate scaling of these spectra, us- 
ing the method detailed in W07, an appropriate model background 
spectrum can then be produced. The SAS tools arfgen and rmfgen 
were used to produce appropriate Auxiliary Response File (ARF) 
and Response Matrix File (RMF) files for the source and residual 
galaxy regions. Finally, the counts recorded in adjacent (raw) spec- 
tral channels were summed to give a minimum of 20 counts per 
spectral bin in the final set of spectra. 



4 OVERVIEW OF THE GALACTIC X-RAY PROPERTIES 

4.1 The contribution of luminous point sources 

By application of the spatial mask described earlier, we are able to 
measure the count rates (net of background) in the soft, medium 
and hard band images associated with both the "bright source re- 
gions" and the "residual emission regions" (hereafter we refer to 
the latter as the "residual galaxy"). Table [4] summarises the results 
for each of the five new galaxies in our sample. Here we have con- 
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Table 4. Contribution of point sources to the total X-ray luminosity. 



Galaxy 


Spillover/Area 
Factors (%) 


Component 


(0.3-1 keV) 


L X (10 39 
(1-2 keV) 


ergs ) 
(2-6 keV) 


(0.3-6 keV) 


NGC300 


12/7 


Bright Sources 


0.19 


0.07 


0.25 


0.51 






Unresolved Sources 


[0.02] 


[0.01] 


[0.03] 








Residual Galaxy 


0.05 






0.05 






Total Measured 


24 


07 


25 


56 


M74 


1 1/24 


Bright Sources 


1.4 


1.3 


3.0 


5.7 






Unresolved Sources 


[0.5] 


[0.4] 


[0.8] 








Residual Galaxy 


1.4 


0.7 




2.1 






lOldl IVlCdaUlCLl 


2.8 


2 


3 


7.8 


NGC3184 


9/22 


Bright Sources 


1.2 


1.4 


3.4 


6.0 






Unresolved Sources 


[0.6] 


[0.5] 


[0.9] 








Residual Galaxy 


2.0 


0.5 




2.5 






Total Measured 


3.2 


1.9 


3.4 


8.5 


M51 


12/24 


Bright Sources 


6.8 


3.8 


6.5 


17.1 






Unresolved Sources 


[0.5] 


[0.4] 


[0.8] 








Residual Galaxy 


6.3 


0.9 




7.2 






Total Measured 


13.1 


4.7 


6.5 


22.3 


M83 


9/30 


Blight Sources 


2.5 


2.0 


3.1 


7.6 






Unresolved Sources 


[0.5] 


[0.4] 


[0.9] 








Residual Galaxy 


3.3 


0.7 




4.0 






Total Measured 


5.8 


2.7 


3.1 


11.6 



verted the measured count rates to equivalent fluxes in the band on 
the basis of a specific spectral model and then transposed from flux 
to luminosity using the distances quoted in Table Q] For the bright 
source regions we assume a spectral model consisting of a power- 
law continuum with photon index r = 1.7 modified by the fore- 
ground absorption in our Galaxy. For the residual galaxy we use 
the model which best fits the actual spectrum of this component, 
as derived in i]5.2| The quoted X-ray luminosities (Lx) have been 
corrected for foreground absorption, that is they are unabsorbed 
values. 

The Lx values reported in Table|4]have been corrected for two 
further effects. The first is the spillover of bright source flux into the 
residual galaxy region (discussed earlier). The second correction 
relates to the fact that the bright source mask obscures a sizeable 
fraction of the galaxy disk. To correct for this, we have scaled up the 
residual galaxy component by a factor derived from interpolating 
the azimuthally-averaged brightness distribution (as represented by 
the radial distributions of the X-ray surface brightness discussed 
in EjGj into the masked regions. The spillover fraction and the area 
correction factor applicable to each galaxy are listed in Table|4] 

The results in Table [4] demonstrate that whereas the set of 
very luminous sources represented by our bright source sample 
provide the dominant contribution to the total galactic X-ray lu- 
minosity above 2 keV, at lower energies, and particularly below 1 
keV, the split between the bright source and residual galaxy con- 
tributions is more balanced, albeit with the exception of NGC 300. 
For the latter galaxy, we were able to set a much lower threshold 
for bright source removal than was possible for the other galaxies 
in the sample (on account of the relatively proximity of NGC 300 
- see Table[TJ. In other words a higher fraction of the integrated X- 
ray flux from point sources was resolved in this galaxy. However, 
we show later (see £j6j> that the low surface brightness inferred for 
the underlying diffuse X-ray emission in NGC 300 is mostly likely 



due to the low SFR density in this Galaxy rather than being simply 
a consequence of a more stringent source rejection threshold. For 
M74, NGC 3184, M51 and M83 a broadly similar source detection 
threshold transposes to luminosity thresholds for source exclusion 
ranging from 2.0 xlO 37 erg s" 1 (for M83) up to 6.0x 10 37 erg s _1 
(for NGC 3184). The likely contribution of unresolved sources to 
the residual galaxy emission in each case is investigated below and 
a correction to a fixed luminosity threshold is applied later in the 
analysis (see ij7}. 

As a further exercise it is possible to estimate the likely con- 
tribution to the residual galaxy emission of an unresolved popu- 
lation of somewhat less luminous sources with spectral charac- 
teristics and distribution similar to the sources represented in the 
bright source list. Here we make use of recent results from Chan- 
dra, which show that the discrete source luminosity function ap- 
propriate to the disk regions of spiral galaxies typical takes the 
form of a power-law with a slope i n the range —0.5 to —0.8 (in 
the cumulative form). For example. iTennant et alj d2001) quote a 
slope of —0.5 for the disk sources in M81 , whereas |Dpane et al.1 
d2004l) measure -0. 6 for NGC 3184 an d IPence et al.1 feOOlh re- 
port -0.8 for MIOl. lSoria & Wul d2003h examine the log N- log 
5 relation for relatively faint HMXBs in the disk of M83 and find 
a broken power-law provides a rea sonable fit to t he d ata, with a 
low-luminosity index of —0.6. Also lColbert et alj d2004l) measure 
a similar range of power-law slopes for the source luminosity func- 
tion within a relatively large sample of spiral galaxies. The flatness 
of this form means that relatively small numbers of very luminous 
sources provide the bulk of the integrated luminosity residing in 
discrete sources. Furthermore, in the context of providing a rough 
estimate of the contribution of unresolved sources to the residual 
galaxy signal, it implies that we can, in principle, extrapolate below 
our source detection threshold to arbitrarily faint levels (assuming 
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of course that the slope of the luminosity function remains constant 
over an appropriately wide range of source luminosity). 

Here we make the assumption that the slope of the source lu- 
minosity function is —0.6 across the set of galaxies and over a wide 
luminosity range. We then use the number of bright sources actu- 
ally observed in the galaxy from the threshold luminosity up to a 
fixed upper cut-off (set here as Lx= 5 x 10 38 erg s _1 - see Table 
0, to define the normalization of the source luminosity function 
for that particular galaxy. We first calculate the integrated emis- 
sion of unresolved sources in the broad 0.3-6 keV band and then, 
on the basis of the bright source spectral model noted earlier, split 
this into soft-, medium- and hard-band contributions. The entries 
in Table|4]under the heading "unresolved sources" summarises the 
results on a galaxy by galaxy basis. In the soft band, which is the 
focus of our attention in describing the X-ray morphology and the 
linkage of X-ray emission to star formation, the potential contri- 
bution of unresolved, but relatively luminous, source populations 
to the measured residual galaxy emission ranges from as low as 
~ 10% in M51 up to ~ 40% in NGC300. By implication the bulk 
of the residual galaxy luminosity in this band must reside in some 
combination of: (i) truly diffuse X-ray emission associated with 
the galaxy disk presumably resulting from the energy input of su- 
pernova and stellar winds and (ii) the integrated emission of pop- 
ulations of relatively soft spectrum, lower-luminosity sources such 
as supernova remnants, cataclysmic variables, active binaries and 
stellar coronal sources, populations which are not well represented 
in the bright-source (i.e., high Lx) samples. There is also a poten- 
tial contribution from a putative extended galactic halo, which is 
often seen in edge-on galaxies; however the face-on aspect of the 
present galaxies suggests that a very extended halo component will 
be difficult to distinguish from the components confined to the disk. 



4.2 X-ray characteristics of the individual galaxies 

Figures[T]&[2]compare the galaxy morphology observed in soft X- 
rays with that measured in the far ultraviolet FUV {X e ff ~ 2267 
A) band by GALEAE, 

although in the case of NGC 3 1 84 we use the 
XMM-Newton Optical Monitor UVW1 (« 2680 A) image recorded 
in the same observation as the X-ray data (see Table[2} to make this 
comparison. In all cases the the UV images were lightly smoothed 
with a circular gaussian mask with a = 4". In a number of the 
galaxies, particularly M74, M83 and M101, the X-ray emission 
traces segments of the inner spiral arms which are delineated as 
very pronounced features in the FUV images. In all cases the cen- 
tral concentration of the soft X-ray emission and its fall-off with in- 
creasing galactocentric radius is quite well matched to the surface 
brightness distribution measured in the FUV channel. This general 
correlation of soft X-ray and FUV light immediately establishes the 
close connection of a significant component of the X-ray emission 
with recent star-formation in these late-type systems. We explore 
the linkage between the X-ray emission and star-formation in more 
detail in Sj6] Here we briefly comment on each of these galaxies, in 
the context of their extended X-ray emitting components. 



4.2.1 NGC 300 

NGC3 00 is an SA(s)d spiral seen at low inclination (i = 
30° ; Ide Vaucouleurs et all 1 1 99 lh at a distance of 2.0 Mpc 



dFreedman et alj2 001). Many large HII regions are visible through- 
out the galaxy, and there is evidence that the galaxy has experi- 
enced many episodes of star formation. However, star formation 
in the central regions appears to be suppressed in recent times, 
with a small population of stars pr esent in the cen tral arcminute 
of the galaxy with ages below 1 Gyr (Davi dgell998h . There is clear 
evidence of spiral structure, although the contrast between "arm" 
and "off-arm" regions in the galaxy is low. In X -rays NGC300 has 
been studied by R OSAT (|Read & Pietsch 1200 lh and mo r e rece ntly 
by XMM-Newton dCarpano et al.ll2005l) . ICarpano et al.l d2005l) re- 



port the presence of diffuse emission in the galactic disk of NGC 
300 with L x < 10 38 erg s" 1 and temperatures of » 0.2 and » 0.8 
keV. In the current paper we determine the aggregate luminosity 
of the 20 brightest point sources within the central 5' of NGC300 
to be Lx~ 5 x 10 38 ergs -1 , in agreement with earlier studies. 
We estimate the X-ray luminosity of the residual emission to be 
~ 5 x 10 37 erg s _1 (0.3-1 keV). The surface brightness of this 
component is extremely low but appears to trace a central bar-like 
structure and to extend out to at least a distance of 2.5 kpc (4' ). 



4.2.2 M74 

M74 (NGC 628) is an SA(s ) c spir al seen almost face-on (i = 
7° ; IShostak & van der Kruitl Il984t) at a distance of 1 1 Mpc 
dGil de Paz et al.ll2007l) . The nuclear region of the galaxy is char- 
acterised by strong HII emission with no evidence for the pres- 
ence of an AGN. A grand-design twin-armed spiral structure 
is o bserved. In X-rays M74 has been stu died by both Chan- 
dra dKrauss et al.ll2005h and XMM-Newton dSoria & Konel 120021: 
ISoria. Pian & Ma zzali 2004), albeit with little attention, to date, on 
the diffuse X-ray component. In the present analysis we confirm the 
presence of an extensive discrete source population in M74 with an 
aggregate luminosity in the 20 brightest point sources amounting to 
5.7 x 10 39 erg s -1 . The underlying residual X-ray emission has a 
luminosity of 2.1 x 10 39 erg s _1 . This component traces the inner 
spiral arm structure and is evident out to a distance of 9 kpc (3' ), 
at which point it falls below the background level. 



4.2.3 NGC 3184 

NGC3184 is an SAB(rs)cd barred spiral seen at low inclination 
(i < 24° ; Lyon-Me udon Extragalactic D atabase [LEDAj3) at a 
distance of 11.6 Mpc dLeonard et al.|[2002h . The nuclear region of 
the galaxy is domina ted by HII emissio n with no compelling ev- 
idence for an AGN dDoane et alj|2004l) . Grand-design twin spi- 
ral arms are observed, although the contrast in brightness be- 
tween "arm" and "off-arm" regions in the galaxy is low. In X-rays 
NGC3184 has been studied by Chan dra with an emphasis on both 
the point source population i Colbert et al.l2004 ) and t he diffuse X- 
ray emitting gas dDoane et alJl2004l ; lTvler et alj|2004l) . The diffuse 
X-ray emissi on shows a stron g correlation with re gions which are 
bright in Ha dTvler et al.l2004l) . lDoane et al.l d2004l) report the pres- 
ence of diffuse X-ray emission concentrated in areas of younger 
stellar populations and star-forming regions, with a surface bright- 
ness five times greater in spiral arm regions than in off-arm re- 
gions. The spectrum of the diffuse emission in the galactic disk 
can be well fitted with a two-temperature th ermal model with kT of 
0.13 keV and 0.43 keV dDoane et al.l2004l) . In the present study we 
measure the aggregate luminosity of the 18 brightest point sources 



3 The GALEX images were obtained from the public archive at 
http://galex.stsci.edu/GR2/ 



4 at http://leda.univ-lyonl.fr 
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to be ~ 6 x 10 ere s , consistent with earlier Chandra estimates 



dColbert et aU2004h . The residual X-ray emission has an X-ray lu- 
minosity of 2.5 x 10 39 erg s _1 and extends from the nucleus out 
to a distance of 8 kpc (2.5' ). 



4.2.4 M51 

M51 (NGC5194) 
on H = 20° 



is an SA(s )bc spiral seen 
Tullvl 1 19741) at a distance 



almost 
of 8.4 



face- 
Mpc 



( Feldmeier. Ciardullo & Jacobvl 1997). It hosts a low-luminosity 



active galactic nucleus, which shows optical emission lines and is 
classi fied as a Seyfert 2 ( Stauffer 1982; Ho . Filippenko & Sargentl 
1997). M51 is interacting tidally with a companion SBO/a 
galaxy (NGC5195), and both systems are seen to con- 
tain starburst activity. There is grand-design twin-arm spi- 
ral struct ure present. In X-ra ys M51 h as been stud i ed by 
Einstein JPalumbo et al] 1 1985h. ROSAT llMarston et alj 1 19951 ; 
lEhle et al.lll995l) , ASCA dTerashima et al.1 fl998h and more re- 
cently by Chand ra jTerashima & Wils on 2004; C olbert et alj2004l ; 
iTvler et alj2004l) and XMM-Newton jDewangan et alj2005l) . Stud- 
ies hav e also been carried ou t at higher X-ray energies by Bep- 
poSAX jFukazawa et"al]|200lh . The diffuse X -ray emission is well 
correlated with the radio continuum and mid-infrared emission dis- 
tributed in the spiral arms establishing a c le ar link with star for- 
mation ( Ehle et al.1 1 19951 ; ITvler et alj |2004|) . iTerashima & Wilson! 
d2004l) find that the diffuse X-ray emission in the southern ex- 
tranuclear cloud is well fitted with a single-temperature thermal 
model with kT of 0.58 keV. The current analysis confirms the 
general properties of the galaxy established by the ROSAT and 
Chandra observations. The aggregate luminosity of the brightest 
25 point sources in M51 is 1.7 x 10 40 erg s" 1 (0.3-6 keV), 28% 
of which is attributable to the nuclear source. Our analysis reveals 
residual soft X-ray emission extending from the bright nucleus 
along the spiral arms out to a distance of 5 kpc (2.5' ). Lower sur- 
face brightness emission from the extended disk is observed out to 
a radial distance of 9 kpc in all directions, with further X-ray emis- 
sion joining M51 to its companion galaxy. The X-ray luminosity of 
the residual component is 7.2 x 10 39 erg s _1 . 



4.2.5 M83 

M83 (NGC 5236) is an SAB(s)c barred spiral seen at low in- 
clination (i = 2 4° ; iTalbot Jensen & Dufourl 1 1979b at a dis- 
tance of 4.5 Mpc IThim et alj 120031) . It harbours an active nu- 
clear region undergoing a violent starburst, a circumnuclear star- 
burst and grand-design twin-armed spiral structure. There are nu- 
merous si tes of active star-f ormation coincident with the spi- 
ral arms JVogler et al] l2005h- In X- rays M83 has been stud- 
ied by Einstein jTrinchieri et al.|[l984l). ROSAT tehle et alJI 19981 ; 
llmmler et all Il999l) . ASCA jOkada. Mitsuda & Dotanil Il997l and 
more recently by Chandra JSoria & Wull2002l ; ISoria & Wull2003l ; 
IColbert et aT]|2004l ; ITvler et alj|2004l) . On an arcminute scale both 
the ROSAT and Chandra X-ray observations trace ap parently dif- 
fuse emission associated with the inner spiral-arms |lmmler et al.1 
ll999l ; ISoria & Wul2003l ; lTvler et alj2004l) . llmmler et"aT] ( ll999h re- 
port the diffuse emission associated with a cluster of bright HII 
regions embedded in the south-western spiral arm as having a two- 
temperature thermal characteristic with kT of 0.26 keV and 0.95 
keV. In the present analysis we measure the aggregate luminosity 
of the bright point sources in M83 to be Lx = 7.6 x 10 39 erg s -1 , 

' 39 erg s _1 from the cen- 



from the bright nuclear region out along the twin spiral arms to a 
galactocentric distance of about 4 kpc (3' ). Extended X-ray emis- 
sion of relatively low surface brightness can then be further traced 
out to ~ 8 kpc. The X-ray luminosity of the residual component is 



4 x 10 39 



erg s 



which includes a contribution of 3.1 x 10 
tral region. The XMM-Newton soft band image traces emission 



5 SPECTRAL ANALYSIS 

Using the methodology outlined in $3.3\ spectral datasets were ob- 
tained for both the bright-source and the residual-galaxy regions for 
four galaxies, namely M74, NGC3184, M51 and M83. For simplic- 
ity, we based our analysis on a single observation for each galaxy, 
namely the observation with the deepest pn exposure (see Table^. 
Due to limited signal-to-noise in the MOS channels, we focus on 
the spectra derived from the EPIC pn camera, except in the case 
of M83, where both pn and MOS spectra were fully utilised and in 
M74 where both datasets were employed in the spectral investiga- 
tion of its bright-source region. Unfortunately the residual-galaxy 
emission in NGC 300 was too faint for detailed spectral analysis to 
be merited. The spectral fitting was carried out using the software 
package XSPEC Version 12.4. 



5.1 Spectra of the bright-source regions 

The spectra representative of the bright-source regions in M74, 
NGC3184, M51 and M83 were well fitted with either a simple 
power-law continuum or a power-law continuum plus an additional 
solar-abundance thermal plasma component (the MEKAL model in 
XSPEC). The best-fit photon indices were all in the range 1.55-1.9 
with the temperature of the thermal emission typically ~ 0.65 keV. 
Fig.|4]shows the measured spectra, the best-fitting model and resid- 
ual x 2 to these best fits. Similarly, Table [5] summarises the details 
of the best-fitting models. After correcting for the flux spillover be- 
yond the spatial mask, the inferred luminosities in the integrated 
bright-source spectra were found to be in general agreement with 
the estimates quoted previously based on the image analysis. The 
photon index observed is typical of the spectra of H MXB and 
LMXB found in the Local Group jGrimm et af]|2007h . The ther- 
mal emission associated with the bright-source regions is of similar 
temperature to the "hot" component observed in the residual-galaxy 
emission (see below). However, where this component is detected 
in the bright-source spectra, namely in NGC3184, M51 and M83, 
it appears at a proportionately higher level than might be inferred 
based on a simple area scaling with respect to the residual-galaxy 
emission. 



5.2 Spectra of the residual-galaxy regions 

Previously we estimated that 9-12 per cent of the counts from the 
bright sources overspill into the residual-galaxy regions. In car- 
rying out the spectral analysis of the latter we have correct for 
this spillover effect, by including an appropriate fraction of the 
best-fit bright-source model as a fixed component of the spectral 
model. Visual inspection of the residual-galaxy spectra revealed 
almost no emission above 1.5 keV in any of the galaxies. It was 
further found that thermal models comprising either one or two 
solar-abundanc^El MEKAL components provided a reasonable fit 



5 In general the quality of the spectra were insufficient to allow useful con- 
traints to be placed on the metal abundances in the thermal plasmas. As 
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Figure 4. The EPIC spectra for four galaxies. Left-hand panels: The spectra of the bright-source regions. Right-hand panels: The spectra of the residual-galaxy 
regions. In all cases the solid line corresponds to the best-fit spectral model (see text). The \ 2 residuals with respect to the best-fitting model are also shown. 
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Table 5. Parameters of the best-fitting models to the spectra of the bright-source and residual-galaxy regions. 
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Galaxy 


Region 


Power-law Cont. 
Photon Index 
Normalization 


Cool MEKAL 

keV 
Normalization 


Hot MEKAL 

keV 
Normalization 


Goodness 
of Fit 
X 2 /dof 


CoolHot 
Flux Ratio 
(0.3-2 keV) 


M74 


Bright Sources 


1.64±0.06 
5.62 x 1CT 5 






146/150 






Residual Galaxy 


- 


0.26±0.04 
3.74 x 10" 5 


[0.65]° 
[< 5.5 x lO" 6 ]" 


197/202 


> 5.4 


NGC3184 


Bright Sources 


1.57±0.10 
4.41 x 10~ 5 


- 


0.71 ±0.1 3 
6.63 x 10" 6 


86/78 


- 




Residual Galaxy 


- 


0.20±0.06 
1.85 x icr 5 


[0.65] a 
[< 3.5 x 10~ e >] a 


122/105 


> 4.1 


M51 


Bright Sources 


1.57±0.05 
2.65 x 1CT 4 


0.19±0.02 
1.20 x 10" 4 


0.61±0.05 

1.48 x 10" 4 


389/378 






Residual Galaxy 


2.90±0.15 
1.14 x icr 4 


0.24±0.02 
1.63 x 10~ 4 


0.64±0.04 
9.94 x 10~ 5 


210/260 


1.3 


M83 


Bright Sources 


1.86±0.02 
5.42 x 10~ 4 




0.61±0.03 
2.78 x 10" 4 


796/679 






Residual Galaxy 


2.92±0.15 
1.12 x 10~ 4 


0.24±0.02 
2.32 x 10" 4 


0.64±0.03 
1.19 x 10" 4 


404/356 


1.5 



" - Upper bound on the normalization of a 0.65-keV thermal component (95% confidence interval). 



Table 6. Physical properties of the diffuse gas present in each galaxy. 



Galaxy Radius" Component Electron Density Thermal Energy Cooling Timescale 



kpc (keV) 10" 3 ??- 1 / 2 cm" 3 lO 55 ^ 1 / 2 erg 10V /2 yr 

M74 10.0 0.2 3.1 2.5 6.6 

NGC3184 10.0 0.2 2.1 1.9 4.3 

M51 7.5 0.2 6.7 2.4 1.9 

0.65 5.2 5.5 8.3 

M83 6.5 0.2 4.9 2.5 4.2 

0.65 3.1 3.6 9.5 



a - Assumed radius of a putative shallow halo component (see text) 



to the residual-galaxy spectra. In all cases a "cool" thermal com- 
ponent was required at kT » 0.2 keV. The spectra with the best 
signal-to-noise ratio, namely those for M51 and M83, also required 
a "hot" thermal component at kT ~ 0.65 keV. These spectral char- 
acteristics appear to be quite typical of th e diffuse components 
seen in normal and starburst galaxies (e.g.,. lFraternali et al„ 2002) 
and consistent with previously publi shed measurements for several 
of the galaxies in the sample (e.g., lEhle et alJI 19981 ; ISoria & Wul 
120031 : iKuntz et alj|2003t Icarpano et al.l 120051 ; W07). The spectra 
of M5 1 and M83 also showed evidence for a "soft excess" which, 
via spectral fitting, was modelled as a soft power-law component 
with photon index ~ 2.9. This latter component contributes 26% 



noted in W07, the fitting of low resolution spectra pertaining to a complex 
multi-temperature plasma with simple models can often result in the artifi- 
cal requirement for strongly subsolar abundances. 



of the flux in the 0.3-2 keV band for M83 and 35% for M51. This 
soft excess may well relate to the contribution of populations of 
soft sources most notably recent supernova remnants and, in fact, 
there is ample observational evidence that such sources are very 
often promin ent in the sp iral ar ms of late-type galaxies, for ex- 
ample IC342 jKong||2003l). M33 jPietsch et alj|2004 [Grimm et al.l 
120051) and M83 dSoria & Wull2003l) . After applying an appropri- 
ate area-correction factor, the luminosities derived from the spectra 
data were broadly consistent with those derived from the imaging 
analysis (see Table|4j. 

There are two points pertinent to the above spectral analysis. 
The first is that the spectral modelling did not require the inclu- 
sion of any absorption over and above the foreground Galactic Nh 
which was included as a fixed component in all the models. The 
implication is that any absorption intrinsic to the host galaxy must 
be relatively low, commensurate with the fact that we are studying 
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systems which have a near face-on aspect. The second point is that 
there was no strong evidence in the residual galaxy spectra for a 
"hard-excess" component attributable to an unresolved population 
of relatively hard discrete sources. This confirms the conclusions 
in £14. 1 1 namely that relatively luminous binaries (albeit below the 
applied luminosity threshold) do not provide a substantial contri- 
bution to of the residual-galaxy emission. In fact the lack of clear 
evidence for hard emission in the spectra of the residual galaxy re- 
gions (and also in the corresponding hard-band images) suggests 
that the estimates in <j4. 1 l of the potential contribution of relatively 
luminous discrete sources might be best considered as upper-limit 
values. 

For each of the galaxies in the sample, we have measured the 
relative contribution of the 0.2-keV and 0.65-keV MEKAL compo- 
nents to the X-ray flux measured in the 0.3-2 keV band. The results 
are summarised in Table [5] where we include the upper limits de- 
rived for a fixed 0.65-keV component in the case of M74 and NGC 
3184. For comparison, the corresponding flux ratio pertaining to 
the cool and hot MEKAL components in M101 (from the spectral 
fitting results in W07) is 3.8. 

If we assume that all of the emission attributed to thermal 
components originates in truly diffuse plasma, then it is possible 
to infer some physical properties of the medium on the basis of a 
particular geometrical configuration. Here we assume that the ther- 
mal plasma is contained within a cylindrical region extending from 
the centre of the galaxy out to 75% of the D25 radius with a half- 
width perpendicular to the plane of 0.5 kpc, representing a shallow 
halo. Using this volume, we can infer the mean electron density n e 
through the derived emission measure TjligV (where r\ is the 'filling 
factor' - the fraction of the total volume V which is occupied by the 
emitting gas). The derived density combined with the volume and 
temperature then leads to an estimate of the thermal energy con- 
tained within the gaseous component. Based on the measurements 
from the spectral fitting, we derive the results summarised in Ta- 
ble [6] In the cases of M51 and M83, where the balance between 
the cool and hot thermal components is most strongly weighted 
towards the hot component, we find that the thermal energy con- 
tained in the hotter component is significantly greater than that in 
the cooler component and that these two components are not in 
pressure balance. Our spectral modelling of the residual emission 
in NGC3184 and M74, although not requiring a 0.65-keV spectral 
component, does not exclude the presence of such a component in 
rough pressure balance with the cooler plasma, consistent with the 
situation pertaining in M101 (W07). The cooling of these plasmas 
is dominated by line emission with radiative timescales in the range 
2 - 10 x 10V /2 yr (Table©. 

Of course, as noted previously, a more likely scenario is that 
the bulk of the observed soft emission originates in the combination 
of diffuse emission closely associated with star-formation in the 
galactic disk and the integrated emission of the populations of soft 
sources within the disk. 



6 THE CONNECTION BETWEEN SOFT X-RAY 
EMISSION AND STAR-FORMATION 

As noted earlier, the general similarity of soft X-ray and FUV im- 
ages establishes a close linkage between the X-ray emission and re- 
cent star-formation in these late-type systems. This connection has 
previously been demonstrated in man y individu a l gala xies and also 
in samples of objects; for example iTvler et all d2004l ) find a very 
strong correlation between X-ray emission and both Ha and mid- 



infrared emission in a set of spiral galaxies observed with Chandra. 
In our early study of M101 (W07), we further investigated how the 
overall point-to-point correlation varies when one contrasts the soft 
X-ray surface brightness with that measured in different wavebands 
extending from the FUV through to optical V band. The fact that 
the best match was between the soft X-ray and U band images was 
interpreted in terms of two spatial components, namely a clumpy 
thin-disk component which serves as a tracer of the spiral arms 
plus a more extended component, possibly located in the lower- 
halo with a larger filling factor, which produces the central concen- 
tration of the soft X-ray emission (W07). However, in the present 
paper we focus on the relationship between the residual galaxy X- 
ray emission and the corresponding FUV emission in terms of their 
radial-averaged spatial distribution. 

Using the masked pn+MOS soft band images (such as in Fig. 
0, we have derived the radial profiles of the soft X-ray emission 
for all six galaxies in our sample. The results are shown in Fig. [5] 
(upper panels). These X-ray radial profiles are based on the sig- 
nal extracted from circular annuli centred on the galactic nucleus, 
with suitable normalization to allow for the regions excluded by the 
spatial masking. Since we wish to compare the X-ray radial distri- 
butions with published FUV and star-formation profiles obtained 
using elliptical annuli (matched to the modest inclinations of these 
galaxies), we have stretched the scale of the X-ray radial axis by 
a factor of typically 5-10%, so as to account for the (marginally) 
different radial gradients obtained if one uses circular as opposed 
to slightly elliptical annuli. 

Figure[5]also shows the corresponding FUV radial brightness 
distrib utions taken from the tabulations in iMunoz-Mateos et all 
d2007h (Table 2 in their paper) except in the case of NGC 3184 
where, since GALEX FUV data are not available, we use XMM- 
Newton Optical Monitor UVW1 (« 2680 A) measurements (with- 
out any internal extinction correction). These GALEX FUV data 
have been c orrected for both foregroun d and internal extinction as 
discussed in Mun oz-Mateos et alj (2007). 

Visual comparison of the soft X-ray and FUV profiles sug- 
gest that radial fall-off seen in the soft X-ray regime is either 
very comparable to that seen in the FUV {e.g. NGC 3184 and 
M83) or somewhat steeper (e .g. NGC 300, M74, M51 an d M101). 
We further use the results of Munoz-Mateos et al. I (f2007h to com- 
pare our current soft X-ray measurements (after applying appro- 
priate scaling factors to convert the X-ray measurements into lu- 
minosity per unit disk area) with the inferred radial dependence 
of the SFR density (i.e. SFR per unit disk area). The results are 
also shown in Fig. [5] (middle panels). We find that in the inner 
disks of our sample of late-type galaxies, the soft X-ray luminos- 
ity generated per unit SFR (hereafter the X-ray/SFR ratio) ranges 
from 2.5 x 10 38 erg s" 1 ( M yr" 1 ) -1 in the case of NGC 300 
up to 4 x 10 38 erg s" 1 ( M yr -1 ) for M74 and finally up to 
~ 10 39 erg s" 1 ( Mq yr- 1 )- 1 forM51, M83 and M101. Note that 
these measurements refer to the soft-band 0.3-1.0 keV luminosity 
after excluding the contribution of the most luminous point sources. 
There is also the suggestion that the X-ray/SFR ratio declines with 
galactocentric radius, at least in the case of NGC 300 and M51. 

Whereas some correction has been applied to the FUV data 
to allow for the extinction intrinsic to the host galaxy, we have not 
applied a similar correction to the soft X-ray data. Earlier we estab- 
lished that based on the spectral fitting, there was no clear evidence 
for any soft X-ray absorption intrinsic to the host galaxies. Any in- 
trinsic soft X-ray extinction would be most pronounced in the inner 
galaxy regions potentially giving rise to a central down-turn in the 
soft X-ray radial distributions - but again there is no compelling 
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Figure 5. A comparison of various radial dependencies for each of the six galaxies in the full sample. In these plots the x-axis refers to the major-axis radius 
scaled to kpc using the assumed distances quoted in Table[T] The vertical dashed line represents the extent of the X-ray extraction region scaled to kpc. For 
each galaxy, the following information is provided: Top pan el: The soft X-ray surface b rightness versus radius (upper curve). The radial profile of the FUV 
emission after correction for internal extinction - taken from Munoz-Mateos et al. 1 2007) (lower curve). Middle panel: The ratio of the soft X-ray luminosity 
in erg s _1 pc — 2 (0.3-1 keV) to the local SFR in units of Mq yr^ 1 pc~ 2 . Bottom panel: Variation in X-ray spectral hardness, (H-S)/(H+S), versus radius, 
where H refers to the 0.8-1.2 keV band and S to the 0.3-0.8 keV band. 



14 R.A.Owen & R.S. Warwick 



evidence for such an effect As a final check we have looked for 
spectral trends versus radius by repeating the imaging analysis for 
two soft X-ray sub-bands, namely 0.3-0.8 keV and 0.8-1.2 keV. 
We then determine a spectral hardness ratio as (H-S)/(H+S), where 
H refers to the 0.8-1.2 keV band and S to the 0.3-0.8 keV band. 
Figure [5] (lower panels) shows the measured variations in spectral 
hardness versus radius. These curves show no marked spectral vari- 
ations towards the inner galaxy regions, thus confirming that intrin- 
sic absorption does not strongly influence the soft X-ray measure- 
ments. 

The demarkation of the above sub-bands at 0.8 keV corre- 
sponds roughly to the position in the M5 1 and M83 spectra at which 
the cool and hot thermal components contribute equally. It follows 
that the spectral hardness profiles in Fig. [5] track the ratio of these 
two thermal components as a function of galactocentric radius. If 
we consider the three sources with best signal-to-noise, we find no 
variation in spectral hardness versus radius in Ml 01, a hint of spec- 
tral softening versus radius in M83 and the hint of an opposite trend 
in M51. 



7 DISCUSSION 

A major objective of our investigation is to intercompare the X-ray 
properties and other characteristics of the extended disks of the six 
late-type galaxies which comprise our sample. To that end we have 
summarised some of the results from the current analysis in Table 
|7] In brief this table provides the following information. We first de- 
fine the extent of the disk region considered in X-rays. For this in- 
tercomparison we have chosen to exclude a central 1' radius region 
in both M51 and M83, so to avoid the complications associated 
with the AGN and/or the nuclear starburst. We have next calculated 
the total SFR within the specified disk region by integrating under 
th e exponential fits to the starb urst density radial profiles reported 
by Munoz-Mate os et al] d2007l) . We next report the luminosity as- 
sociated with the disk by taking our soft band residual luminos- 
ity estimates and then correcting both to a common source exclu- 
sion threshold of Lx> 10 37 erg s _1 (0.3-6 keV) (using essentially 
the same methodology as applied in S|4.2t and to a broader 0.3-2 
keV band-pass (using the spectral models for the residual emission 
regions discussed earliefl). These corrections were generally less 
than 20% except for NGC300, where the luminosity increased by a 
factor of 2.2 as a result of a significant additional contribution from 
discrete sources. Finally, we also note the relative strengths of the 
cool (0.2-keV) and hot (0.65-keV) thermal components in terms of 
the ratio of their fluxes in the 0.3-2 keV (band), as inferred from 
the spectral modelling ( i]5.2b . 

6 We can use t he F UV extinction estimates employed by 
iMunoz- Mateos et al] fe007t) to predict the soft X-ray absorption us- 
ing the scaling A(FUV)/N H = 1.6 X 10" 21 mag cm" 2 , which 
assumes a Galactic dust- t o-gas r atio and extinction law with Ry = 3.1 
iBohlin. Savage & Drake! 1 19781 ; ICardelli, Clayton & Mathisl 1 19891) . This 
leads to estimates of the intrinsic X-ray absorption in the inner disk regions 
in the range Nh = 5 — 12 X 10 20 cm -2 , across the set of galaxies. 
Since these values are too high to be consistent with spectral fitting we 
can conclude either that the assumptions underlying the FUV to X-ray 
extinction scaling are not valid or, as seems likely, the FUV sources 
are more strongly embedded in obscuring material than the soft X-ray 
emission. 

7 In the case of NGC300 we assume a single component (0.2-keV) 
MEKAL spectrum and for M101 we use the spectral model reported in 
W07. 



The results summarised in Table [7] establish that averaged 
over the defined disk regions, the X-ray/SFR ratio varies from a 
maximum of ~ 10 39 erg s _1 ( Mq yr" 1 )" 1 in the case of M51 
and M83 down to w 2 x 10 38 erg s" 1 ( M Q yr" 1 )" 1 in M74. 
However, these disk-average values mask a trend, which is evident 
in Fig. [6] where we have plotted the X-ray/SFR ratio versus the 
local SFR density. Above a local SFR density of ~ 3 x 10~ 8 
Mq yr -1 pc~ (corresponding to the inner disks of M51 and 
M83), the X-ray/SFR ratio is consistent with an upper bound value 
of 1 — 1.5 x 10 39 erg s -1 ( M yr -1 ) -1 . However, below a lo- 
cal SFR density of ~ 3 x 10~ 8 M yr" 1 pc" 2 , there is signif- 
icant scatter in the ratio, spanning a range up to a factor 5 below 
the maximum. On the one hand some of this scatter appears to be 
due to the fact that all five galaxies show evidence for a fall-off 
in the X-ray/SFR ratio in their outer disks. On the other hand sys- 
tematic errors in both the X-ray and FUV measurements will grow 
in importance as one tracks a signal of rapidly decreasing surface 
brightness against an essentially constant background. Neverthe- 
less, these results do suggest that the efficiency of X-ray production 
may be sensitive to the local level of star-formation activity in the 
disk. 

In a recent study, Ma s-Hesse, Oti-Floranes & Cervind J200g) 
(hereafter MH08) have investigated the temporal evolution of both 
the soft X-ray and far infrared (FIR) luminosity following a burst 
of star formation. The starting point of their analysis is the adoption 
of a stellar population synthesis model which follows the evolution 
of a cluster of massive stars formed either at the same time (which 
MH08 refer to as an instanteous burst, IB) or at a constant rate over 
a period of several tens of Myr (an extended burst, EB). In either 
scenario, diffuse soft X-ray emission is produced through the heat- 
ing of bubbles within the interstellar medium to temperatures of 
~ 10 6 " 7 K as a result of the mechanical energy input from the 
winds of massive stars created in the starburst or the eventual de- 
struction of such stars in supernovae. Their model also includes the 
soft X-rays emitted by individual SNR during the adiabatic expan- 
sion phase of their evolution, but excludes any contribution from 
the stellar atmospheres (which in any case would be extremely 
small) or accreting high-mass X-ray binaries created in the region. 
MH08 (see their Figure 1) find that the soft X-ray luminosity in- 
creases rapidly during the first few Myr following the onset of star- 
burst activity, largely through the input power of the stellar winds 
of the most massive stars. After the first 3 Myr, the most massive 
stars produce supernovae and thereafter it is the mechanical energy 
of such explosions which dominates the energy budget. In the IB 
models the soft X-ray luminosity peaks at ~ 3 Myr, whereas in the 
EB models it continues to rise, albeit more gradually, until even- 
tually an equilibrium is established between the formation and de- 
struction of stars (after typically ~ 40 Myr for solar metallicities). 

In comparing our present measurements with the predictions 
of MH08 it is important to bear in min d that we have used SFR es - 
timates based on FUV measurements ( Mun oz-Mateos et alj|2007t) . 
which in t urn re ly on the UV-to-SFR cal ibration su ggested by 
iKennicuttl dl998h . As noted by MH08, the iKennicuttl j |l998t) cal- 
ibration assumes star formation over a wide 0.1-100 Mq mass 
range, whereas the SFR values they quote are for a more restricted 
2-120 Mq range (which has a more direct bearing on the result- 
ing soft X-ray and FIR luminosity). Assuming a Salpeter initial 
mass function, one obtains the scaling SFR(2-120 Mq)/ SFR(0.1- 
100 Mq) = 0.293, implying that our estimates of the X-ray/SFR 
ratio should be increased by a factor 3.4 when making a comparison 
with the MH08 predictions. Our upper-bound estimate, after apply- 
ing the above scaling, becomes ~ 5 x 10 39 erg s _1 ( Mq yr -1 )" 1 , 
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Table 7. Derived properties of the galactic disks. 



Galaxy 


Disk Radial Extent 
Inner (kpc) Outer (kpc) 


SFR 

M yr- 1 


L x (0.3-2keV) a 
10 39 ergs" 1 


CookHot 
Flux Ratio 


NGC300 





2.9 


0.19 


0.11 




M74 





16.7 


5.2 


1.2 


> 5.4 


NGC3184 





12.2 




3.4 


> 4.1 


M51 


2.4 


8.9 


7.5 


6.8 


1.3 


M83 


1.3 


8.5 


3.5 


3.4 


1.5 


M101 





20.9 


8.1 


4.3 


3.8 



a - Luminosity after excluding sources with Lx> 10 37 erg s 1 (0.3-6 keV) 




loq(SFR density) 

Figure 6. The ratio of the soft X-ray luminosity in units of erg s — 1 pc~ 2 (0.3-2 keV) to the local SFR density in units of Mq yr — 1 pc — 2 versus the local 
SFR density. The data points correspond to NGC300 (orange squares), M74 (green circles), M51 (red triangles), M83 (black points) and M101 (blue crosses). 



which correponds to the X-ray/SFR ratio predicted some 10 Myr 
after the onset of an extended burst of star-formation, assuming a 
1% efficiency in the conversion of mechanical energy into X-ray 
emission. Although it can be argued that star formation in the disks 
of late-type spiral galaxies generally proceeds over long-periods of 
time, the underlying activity pattern is presumably one of localised 
bursts of star formation, so a match with "young" extended bursts 
is perhaps not unreasonable. Of course, in matching the observa- 
tions to the predictions, the efficiency of X-ray production and the 
duration of the EB episode are inversely linked, for example the 
above X-ray/SFR ratio is also reached after 30 Myr in the EB sce- 
nario, provided the energy conversion efficiency is scaled down to 
just 0.2%. Clearly a more detailed study, taking into account addi- 
tional information relating to the current star formation rate and the 
star formation history relevant to each galaxy disk, might well pin 



down these parameters - but such an analysis is beyond the scope 
of the current paper. 

A number of recent studies have investigated whether the 
integrated soft (nominally 0.5-2 keV) or hard (2-10 keV) X-ray 
luminosities measured in galaxies not dominated by a central 
AG N might serve as a proxy for the underlying galaxy-wide SFR 
(e.g. Grimm, Gilfanov & Sunyaev 2003; Ranalli, Comastri & Setti 
l2003LlGilfanov. G rimm & Sunyaev 20 o3), which might then be ap- 
plicable to mor e distant objects {e. g. Rosa-Gonzal ez et al .120071) . If 
we consider the lRanalli. Comastri & Settil j2003h calibration of the 
SFR in terms of the integrated soft X-ray luminosity (their equation 
14), apply the corrections noted by MH08 and also scale the X-ray 
luminosity from a 0.5-2 keV bandpass to the 0.3-2 keV bandpass 
(which increases the luminosity by a factor of 1.1 — 1.5 depend- 
ing on the spectral assumption), we obtain an X-ray/SFR ratio in 
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the ranee 10 - 14 x 10 39 erg s" 1 ( M yr" 1 )" 1 which is fac- 
tor 2 — 2.8 higher than the (upper-bound) value derived from the 
current analysis. The difference is, of course, due to the fact that 
Ranall iTComastri & Settil J2003I) use the integrated soft X-ray lu- 
minosity of the entire galaxy, which includes a very substantial 
contribution from the bright discrete source population. In fact the 
ratio of the total galaxy luminosity to the residual galaxy luminos- 
ity over the 0.3-2 keV band for our current set of galaxies (see Table 
[4]l is typically in the range 2 — 2.5, consistent with th i s pictu re. It 
follows that when using the lRanalli. Comastri & Settil J2003h cali- 
bration to convert soft X-ray measurements to SFR rate estimates, 
an implicit assumption is that both the rate of production of X-ray 
binaries following a starburst and the subsequent balance between 
point source and diffuse luminosity is comparable to that pertaining 
in local galaxies. 

As discussed above there is evidence that the X-ray/SFR ratio 
falls as the local SFR density declines, implying that the interstel- 
lar environment influences either the starburst and/or the efficiency 
with which its energy is reprocessed into the X-rays. Factors which 
may potentially be relevant include the local gas density (relevant 
to the conversion of mechanical energy into heat via shocks) and 
the local metallicity (which may, for example, influence the evo- 
lution timescales of massive stars, their stellar wind properties and 
the emissivity of the X-ray plasma heated by the starburst). It is 
also notable that the galaxies in our current sample with the high- 
est X-ray/SFR ratios are those with the strongest weighting towards 
hot plasma (as measured by the coohhot flux ratios in Table|7J. In 
this context it would seem entirely plausible that X-ray production 
efficiency and "mean" plasma temperature are coupled quantities. 
Although the current data are too restricted to merit further con- 
sideration of these issues, such effects may well be relevant when 
attempting to fine tune the calibration of the soft X-ray luminos- 
ity versus SFR relation for specific applications in which the bright 
source component is spatially resolvable. 

Star formation in the disks of spiral galaxies is thought to 
be triggered by the passage of a spiral density wave through 
the ISM. The resultant massive star formation results in copious 
FUV emission which in turn is rapidly reprocessed into Ha and 
FIR emission. As discussed above, the mechanical energy input 
from stellar winds and supernovae eventually gives rise to X- 
ray emission. This would seem to give an immediate explana- 
tion of the strong correlation observed between the infrared, op- 
tical and UV star-formation t racers and the soft X- r ay emission 
seen in spiral galaxie s (e.g., iRead & Pietschl 1200 ll ; iKuntz et"al] 
l2003l ; lTvleretalj2004 W07). However, the underlying timescales 
are not identical. UV and Ha emission will be produced on a 
timescale commensurate with the lifetime of the most massive 
stars (~ 3 x 10 6 years), whilst according to model predictions 
dLeitherer & Heckman|[l995l ; ICervino. Mas-Hesse & Kunthll2002l ; 
MH08), it takes approximately ten times longer to maximise the 
diffuse X-ray emission. Given the delay between the peak in 
UV/Ha emission and the X-ray heating, one might predict a spatial 
offset between the diffuse X-ray emission and the other spiral trac- 
ers more closely tied to the passage of the spiral density wave. For 
example, W07 argue that in the case of M 101 one might predict a 
rotational lag of » 24° at a galactocentric radius of r ~ 7.5 kpc 
(given a rotational velocity v rot ~ 200 km s" 1 and assuming that 
the spiral pattern corotates with the disk material at r ~ 15 kpc). 
However, lags of this magnitude are n ot generally evide nt in the 
X-ray morphology of late-type galaxies dTvler et al.l2 004). 

Assuming a gas filling factor 77 ~ 1 and a spatial extent per- 
pendicular to the plane of the galaxy of 0.5 kpc, we earlier derived 



radiative cooling timescales in the range 2 — 10 x 10 yr (Table 
[6). These estimates are comparable to the galaxy rotational periods 
in the region sampled by the current X-ray measurements and are 
clearly inconsistent with the presence of narrow X-ray spiral arm 
features. Matching the radiative timescales to the "young" extended 
burst scenario discussed earlier would require 77 <<1, consistent 
with a clumpy thin-disk distribution for the X-ray emitting plasma 
most closely associated with the spiral arms. In practice, however, 
the observed X-ray emission will encompass a complex web of dif- 
fuse features extending over a range of spatial scales including su- 
perbubbles and, i n the most active region s , outflows into the lowe r 
galactic halo (e.g. IStrickland et al.ll2004al : IStrickland et al.ll2004bb . 
Where there is no confinement by chimneys or similar structures 
in the ISM, energy losses arising from adiabatic expansion of the 
hot gas in the disk may also help localise the spiral arm compo- 
nent (W07). Relatively luminous discrete sources including young 
SNRs, some perhaps linked to disk population stars in interarm re- 
gions, together with aggregations of less luminous objects, add fur- 
ther to the mix. 



8 CONCLUSIONS 

We have used XMM-Newton observations in a study of the extended 
extranuclear disks of six nearby face-on spiral galaxies. Using a 
novel spatially masking technique to minimise the impact of the 
brightest discrete sources, we have investigated both the spatial 
morphology and the spectral properties of the residual soft X-ray 
emission emanating from the disk regions. The strong correlation 
found between soft X-ray and FUV images in terms of the tracing 
of specific spiral features and the overall radial extent of the emis- 
sion unambigously establishes a close link between the soft X-ray 
emission and recent star formation. 

More detailed comparison of the radial profiles of soft X- 
ray and FUV surface brightness distributions establishes an X- 
ray/SFR ratio of ~ 5 x 10 39 erg s" 1 ( M Q yr" 1 )" 1 (referred to 
the 0.3-2 keV band and where the SFR applies to stars in the mass 
range 2-120 M Q ) for the inner disks of M51 and M83. This is 
roughly a factor 2.5 below the soft X-r ay to SFR calibration re- 
port by iRanalli. Comastri & Settil d2003l) . consistent with the fact 
that our estimate excludes the contribution of discrete sources with 
L x > 10 37 erg s" 1 (0.3-6 keV). Our measured X-ray/SFR ratio 
matches that predicted some 10 Myr after the onset of an extended 
burst of star formation in the models of MH07, when the efficiency 
of the conversion of mechanical energy input to X-rays is set at 
~ 1%. Our observations suggest a fall-off in the X-ray/SFR ratio 
as the local SFR density declines, implying that the interstellar en- 
vironment influences either the starburst and/or the efficiency with 
which its energy is reprocessed into the X-rays. 

With the bright sources excluded, the spectra of the residual 
galaxy regions are, in the main, well matched by a two-temperature 
thermal plasma model. The characteristic temperatures of ~ 0.2 
and ~ 0.65 keV are in line with published results for other spiral 
and starburst galaxies. The relative strengths of these two compo- 
nents varies across the sample, with the galaxies having the highest 
X-ray/SFR ratio characterised by a higher "mean" temperature. The 
physical properties of the gas found in the disk are shown to be con- 
sistent with a clumpy thin-disk distribution presumably composed 
of diffuse structures such as superbubbles together with individual 
SNRs and other source aggregations. 

Future observations should allow more detailed investigation 
of a wide range of factors, for example the local gas density and 



X-ray emission from spiral galaxies 17 



the local metallicity, which might influence the X-ray/SFR ratio 
measured in late-type galaxy disks. 
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